Determination of the antiferroquadrupolar order parameters in UPds 
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By combining accurate heat capacity and X-ray resonant scattering results we have resolved 
the long standing question regarding the nature of the quadrupolar ordered phases in UPda. The 
order parameter of the highest temperature quadrupolar phase has been uniquely determined to 
be antiphase Qzx in contrast to the previous conjecture of Q^Q._yi. The azimuthal dependence of 
the X-ray scattering intensity from the quadrupolar superlattice reflections indicates that the lower 
temperature phases are described by a superposition of order parameters. The heat capacity features 
associated with each of the phase transitions characterise their order, which imposes restrictions on 
the matrix elements of the quadrupolar operators. 

PACS numbers: 75.25.-(-z, 75.10.-b, 75.40.Cx, 78.70.Ck 



In the past couple of decades there has been consider- 
able interest in the orbital orderin g of d and / electron 
systems HiSSiiHiSMIIll The highly 
degenerate /-electron shells in actinide and rare earth 
systems provide a wealth of local degrees of freedom: 
dipolar (magnetic), quadrupolar, octupolar etc. In lo- 
calised /-electron systems these degrees of freedom be- 
come potential order parameters, which can lead to in- 
teresting and complex phase diagrams. Whilst, histor- 
ically, the order associated with magnetic dipole mo- 
ments has been studied extensively, more recently the 
importance of electric quadrupoles in magnetic materi- 
als has been recognized |12i . In classical electrodynamics 
the multipole expansion suggests that the interaction be- 
tween higher order multipoles is seemingly much weaker 
than between dipoles. However, the interaction between 
multipoles is quantum mechanical in origin, and dipo- 
lar and quadrupolar interactions may be equally strong. 
More recently, quadrupolar order has been associated 
with new and interesting behaviour, such as the novel 
heavy Fermion state in PrFe4Pi2 6], and the exotic su- 
perconductivity in PrOs4Sbi2i which may be mediated 
by quadrupolar fluctuations jJI- 

UPds is a very interesting system, being a rare example 
of a localised uranium intermetallic compound, as well as 
belonging to the small class of metallic materials which 
exhibit long-range quadrupolar order. Despite intensive 
experimental investigation over the past 25 years an un- 
derstanding of the series of four phase transitions below 
8 K, and the exact nature of the quadrupolar ordering 
in UPda have proved highly challenging. For the first 
time we are now able to distinguish which order parame- 
ter is associated with the highest temperature quadrupo- 
lar phase between T — 7.8 K and 6.9 K. This has been 
achieved using the unique properties of X-ray resonant 
scattering (XRS), which may be visualised as a process in 



which an incident photon promotes a core electron to an 
excited intermediate state that then decays back to the 
initial state, emitting a scattered photon. XRS reveals in- 
formation about the ordering of the ground state of mul- 
tipoles associated with this intermediate state, since at 
resonance the scattering length becomes a tensor Q| 
whose elements are directly related to the multipole mo- 
ments. Quadrupolar order describes the periodicity of 
charge distribution asphericities and additional super- 
lattice reflections are observed in the case of antiferro- 
quadrupolar (AFQ) ordering. Characterisation of the 
dependence of such reflections on photon energy and po- 
larisation, and azimuthal angle (rotation around the scat- 
tering vector) allows details of the multipolar structure 
to be deduced by comparison with calculations based on 
tensors for the different order parameters 

UPds crystallizes in the double-hexagonal close-packed 
(dhcp) structure (lattice parameters a = 5.73 A, c = 9.66 
A), with uranium ions sitting at sites with locally hexag- 
onal and quasi-cubic symmetry (space group P63/mmc). 
Four transitions at Tq = 7.8 K, T+i = 6.9 K, T_i = 6.7 
K and Ti = 4.4 K have been revealed by a combina- 
tion of microscopic e.g. neutron [la llfl Il2l and X- 
ray scattering, and macroscopic measurements: ul- 
trasound 1 191 . heat capacity ^20, 21], magnetic suscepti- 
bility [21I |22|. thermal expansion and magnetostriction 
|23| . The quadrupolar phase transitions are very sensi- 
tive to doping Np for 111241 , and Pt for Pd 23] . Neutron 
MfilhTlhd and X -ray [ll| scattering studies indicate that 
the phase transition at Tq is to an AFQ structure asso- 
ciated with periodic lattice distortions and a doubling 
of the double-hexagonal unit cell leading to superlattice 
reflections at Q = (2/i -|- 1, 0, 1) in the orthorhombic no- 
tation, which will be used hereafter. Polarised neutron 
diffraction (PND) measurements 0] suggested that the 
phase between Tq and T+i was Q^2_y2, and the earlier 



2 



X-ray experiments jllf were consistent with this hypothe- 
sis. However, due to Hmitations on cryostat technology, it 
was not possible to perform X-ray azimuthal scans at that 
time in the appropriate temperature regimes. Recently, a 
new crystal field model |25| with a doublet ground state 
on the quasi-cubic sites, as opposed to the earlier sin- 
glet ground state model of Buyers et al. , provided a 
qualitative understanding of the succession of quadrupo- 
lar phase transitions with suggestions for the possible 
order parameters. It also presents a framework for un- 
derstanding whether the transitions are first or second 
order. In this paper, we present X-ray resonant scat- 
tering data which shows unequivocably that below Tq 
the AFQ structure is described by the Qzx, rather than 
Qx2-y2, order parameter. 

Previous heat capacity experiments were made before 
the T+i and T_i transitions had been identified, and 
so we have undertaken new measurements to distinguish 
the two transitions. The heat capacity of polycrystalline 
samples of 51.33 mg UPda and 16.79 mg ThPda, used 
as a phonon blank, was measured from T — 2 — 300 
K using a PPMS-9 Quantum Design calorimeter in the 
Actinide UserLab at the Institute for Transuranium Ele- 
ments. Contributions to the heat capacity from the sam- 
ple holder and grease were measured separately and then 
subtracted from the total signal. 

Heat capacity results reveal a lambda anomaly below 
r_i, which is clearly first order, whilst there are no large 
entropy changes at the other transitions, see Fig.^ This 
affects the matrix elements for the order parameters. Us- 
ing the three-level model developed by McEwen et al. 
[25] , the matrices representing the quadrupolar operators 
can be written as 

/ A A\ / A' A' \ 

Q,2_y2 ^ { A B ]q,^ = { A' B' ] (1) 
\AB0j \A'B'oJ 
/ -Ai At \ / A'i -A'i \ 
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where the A^'^ terms mix the singlet with the doublet 
states, and the -B^'-' terms split the doublet. From Lan- 
dau theory the operators' symmetries mean that there 
should be at least two first order transitions, while the 
data in Fig. ^ indicates that there is only one strongly 
first order transition at r_i. Since the B^'^ term splits 
the ground doublet, leading to entropy changes at the 
transition, either B or B' must be ~ 0, to make Q = ~Q 
for either QxS-yS or Qzx, such that there can be just one 
strongly first order transition. 

For the XRS experiment a single crystal of UPda 
was grown, at the University of Birmingham, using the 
Czochralski method with starting materials of 3N U and 
4N Pd. The experiment was performed on the ID20 
beamline at the ESRF. Resonant scattering studies can 
be carried out at the L and M edges of uranium com- 
pounds, but since we wish to probe the / electrons, we 
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FIG. 1: a) Heat capacity of UPda and ThPda (open blue 
squares and a black line respectively), up to 200 K at which 
point the phonons are the dominant contribution to the heat 
capacity, b) UPda heat capacity in the region of the four 
transitions. The lambda anomaly clearly shows that the first 
order transition is associated with T_i rather than T+i. The 
inset shows ultrasonics data Ti] from which the T-i and T+i 
transitions were first identified, with hysteresis in 044 at T_i. 



used the Miy-edge, E = 3.726 keV, at which dipolar 
transitions connect the core 3c?3/2 states to the 5/ states. 
The electric dipole transitions, El, dominate the reso- 
nant scattering cross section. 

At E = 3.726 keV only a small region of reciprocal 
space is accessible. Therefore our crystal was cut with a 
reciprocal space (207) face such that both the (103) and 
(104) superlattice peaks could easily be measured over a 
wide range of ^E". The 417 mg sample was polished with 
0.25 /im diamond paste and mounted in an azimuth dis- 
plex cryostat, allowing us to operate the diffractometer 
in the vertical plane, in which the incident X-rays are 
cr-polarised: see Figure [3 for a schematic of the experi- 
mental set-up. This makes azimuthal rotation possible, 
unlike in the previous horizontal geometry experiment 
|llj | . Polarisation analysis was performed by mounting an 
Au (111) analyser crystal. Data was normalised against 
the monitor. 

The azimuthal dependence of the scattering intensity 
at the (103) refiection at T = 7.1 K, i.e. within the 
first quadrupolar phase, is shown in Figure 01 The az- 
imuth angle 5" is defined relative to the reference vector 
[010]. At each azimuthal angle, rocking curves were taken 
and the intensities of the aa' and (ttt' polarization chan- 



FIG. 2: A schematic of the set up for measuring the azimuthal 
dependence of the scattering intensity of a superlattice reflec- 
tion at Q. cr-polarised X-rays incident at an angle 6 on to the 
sample are scattered by a vector Q. The scattered X-rays are 
separated into the rotated, tt', and unrotated, a' , channels by 
a polarisation analyser before detection. The sample stage 
is rotated to give an azimuthal angle 'if about the scattering 
vector. The unit vectors Ui define the reference frame. 

nels determined by fitting a Lorentzian squared lineshape 
minimizing chi-squared. 

The azimuthal dependence of the allowed order param- 
eters was calculated by summing the second rank tensors 
Tn of the individual uranium quadrupoles to construct 
the resonant scattering length of the unit cell dl] 

/ = ^T„exp(iQ-r). (2) 

n 

The scattering amplitude is then given by 

A = e'-f-e (3) 

where the incident (e) and scattered (e') polarisations 
are transformed into the coordinate system of Blume and 
Gibbs following the method of Wilkins et al. 

FigureElshows that the (103) data is in excellent agree- 
ment with the above calculation for the azimuthal de- 
pendence of Qzx antiferroquadrupolar order. In the air' 
channel, note the asymmetry about = 0° in both the 
data and Qzx calculation, and that the maxima are not 
at \E' = 0,180°. In the aa' channel, the data and cal- 
culation show a broad minimum at 90° with symmetry 
about \1/ = 0°. Since our sample was cut with a (207) 
reciprocal lattice face, the scattering vector Q = (103) 
is not coUinear with the face normal (n). The a' po- 
larisation vector is perpendicular to both n and Q so 
the non-coUinear nature is not observed, resulting in the 
symmetry about 5* = 0° seen in aa' azimuthal measure- 
ments. However, the tt' polarisation vector lies in the 
plane of n and Q, and the non-coUinearity leads to an 
asymmetry in the azimuthal variation in intensity. Calcu- 
lations of the dependence of the scattering intensity for 
the Qx2^y2, Qxy and Qyz order parameters do not agree 



FIG. 3: The azimuthal dependence of a) the crvr' and b) 
the aa' scattering intensities of the (103) peak in UPds at 
the U Miv edge, at T = 7.1 K. Comparison is made with 
calculations for the Qzx (solid line) and the Q^2_y2 (dashed 
line) order parameters. 




FIG. 4: The Qzx AFQ structure with antiphase stacking 
along the z-axis in UPda at T+i < T < To in an orthorhombic 
unit cell. The U 5/ quadrupoles on the quasi-cubic sites are 
represented schematically by ellipsoids. 



with the data, as they show either the wrong periodicity, 
symmetry or maxima and minima positions, and hence 
these order parameters can be ruled out. Qzx provides a 
natural explanation for the macroscopic distortion to the 
orthorhombic cell 23] due to the splitting of the x — y 
symmetry, see Fig. 01 Combining the knowledge that the 
transition at Tq — 7.8 K is to a Q^^ AFQ ordering of the 
5/^ electrons, with the heat capacity evidence that this 
transition is either second order or very weakly first order 
requires B' ~ 0. Although the discussion in '25] assumed 
the order parameter was Qx2_y2, with _B ~ 0, this model 
is equally valid for Q^x as the operator matrices have the 
same symmetry, see equation 

In order to investigate the change in quadrupolar or- 
der in the different phases of UPds seen in various mea- 
surements and the temperature dependence of the RXS 
|lll |. azimuthal measurements were made about both 
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Q = (103) and (104) at T = 5.2 K to study the third 
quadrupolar phase, see Fig.O These measurements were 
made in the cttt' channel only, since there is considerable 
charge scattering interference in the ao' channel. From 
the first order lambda peak in the heat capacity at T_i, 
we expect to find evidence of an admixture of Q^i _y2 in 
this phase. It is very difficult to stabilize the cryostat 
reliably between 6.7 and 6.9 K for extended periods and 
so we have concentrated on measuring azimuthal scans 
in the third quadrupolar phase. However, the heat ca- 
pacity data indicates that at T+i, AS* 0, which would 
be consistent with the evolution of Qyz order. 
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FIG. 5: The azimuthal dependence of the crvr' scattered in- 
tensity from a) the (103) and b) the (104) superlattice peaks, 
in UPds at the U Miv edge, at T = 5.2 K. Dashed lines show 
least squares fits to the data as described in the text. 

The azimuthal dependence of scattering from the (103) 
peak reflects the part of the AFQ structure which varies 
in antiphase along the c-axis, whilst that for the (104) 
peak reflects that which is uniform along the c-axis (in- 
phase stacking) . Clearly the azimuthal dependence of the 
Q = (103) scattering at 5.2 K is more complicated than 
in the higher temperature phase. A least squares fit 
to the data, varying the contributions from the allowed 
quadrupolar moments, indicates that Qzx still predomi- 
nates in the order parameter, but that Qxy and Q-^2_y2 
are also present. The fit in Figure Et) is derived from 
an anomalous scattering tensor in which Qzx, Qxy and 
Q^2_y2 are present in the ratio 77:10:13. The (104) re- 
flection, which is observed below T+i, indicates the ex- 
istence of in-phase (along the c-axis) components to the 
order parameter, as well as the antiphase components de- 
duced from the (103) reflection. The fit to our data in 
Figure ^p), which reproduces the periodicity of the (104) 
data, indicates that the most significant components are 
Qxy and in the ratio 2:1. 

In conclusion, our new heat capacity measurements re- 
veal only one strong first-order transition at T = 
whilst AS ~ at both T+i and Tq. These observa- 
tions constrain the quadrupolar operator matrix elements 



within our model [25j. Our azimuthal XRS experiments 
have unambiguously identified the order parameter for 
the first antiferroquadrupolar phase, T+i < T < To, in 
UPds as Qzx- This may be reconciled with the PND 
(5x2 -y2 result, since those measurements were necessar- 
ily made in a magnetic field. Calculations show that the 
wavefunctions and energies of the crystal field states in 
our model j25l | are significantly modified in fields of a few 
Tesla. We will examine with XRS the field behaviour of 
the different quadrupolar phases in future experiments. 
Our results have provided further valuable insight into 
the nature of the other quadrupolar phases, indicating a 
surprisingly complex sequence of order parameters lead- 
ing to a rotation of the charge densities. 
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